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1.  INTRODUCTION 


An  integral  part  of  the  ID -2D  synthesis  technique  used  to  esti¬ 
mate  the  energy  coupled  into  the  ground  from  a  near -surface  burst  is  a 
simplified  air  fireball  model.  The  current  model  used  in  ground  coupling 
work  was  developed  by  Pyatt,  et  al.  in  1969.  One  of  the  main  limitations 
of  this  model  is  the  assumption  that  the  device  output  is  a  square  wave  in 
time;  i.  e. ,  the  device  radiates  at  a  constant  rate  for  a  finite  period  of  time. 
Analytic  modeling  work  in  the  coupling  problem  has  indicated  that  the  time 
dependence  of  the  device  output  is  an  important  parameter  in  determining 
the  radiative  energy  coupled  into  the  ground.  This  indication,  together  with 
the  fact  that  a  square  wave  is  a  poor  representation  of  the  output  of  a  device, 
led  us  to  reformulate  the  Pyatt  air  model  for  a  more  realistic  device  radia¬ 
tive  yield.  During  the  course  of  this  reformulation,  several  other  improve¬ 
ments  were  made  to  the  earlier  model. 

The  purpose  of  this  report  is  to  describe  in  some  detail  our  re¬ 
formulation  of  the  air  fireball  model.  We  shall  assume  the  reader  is  fa¬ 
miliar  with  the  Pyatt  model  in  general  terms.  As  in  this  early  model,  our 
reformulation  envisions  the  air  "burning  out"  at  a  fixed  temperature.  That 
is,  for  air  temperatures  less  than  the  burnout  temperature,  the  air  opacity 
is  assumed  very  large,  and  for  air  temperatures  greater  than  the  burnout 
temperature  the  air  opacity  is  assumed  very  small. 
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2.  REPRESENTATION  OF  THE  DEVICE  OUTPUT 

We  assume  a  point  source  of  radiation  of  yield  Y  and  a  time 
dependence,  f(t)  ,  given  by 

f /a  _  1  (a+/3)e^t-to^ 

Here  N  is  a  normalization  constant  such  that 

CC 

J  dtf(t)  =  1  . 


This  function  has  the  properties : 

a)  Rises  like  eat  for  small  t  ; 

b)  Falls  like  e"^  for  large  t  ; 

c)  Has  a  single  maximum  at  t  =  t  : 

o 

d)  f(tQ)  =  1/N  . 


3.  BuRNOUT  IN  INFINITE  AIR 


We  assume  that  during  the  burnout  phase  of  the  fireball  growth 
the  air  heats  to  a  temperature  T0  (the  burnout  temperature)  and  gets  no 
hotter.  We  designate  by  Q  the  energy  density  required  to  heat  the  air  to 
this  temperature  Tfi  .  At  this  temperature,  we  assume  the  air  has  a 
small  residual  opacity  x  .  The  fact  that  x  t  0  means  that  a  thermal 
radiation  field  will  build  up.  Let  R  =  R(t)  by  the  fireball  radius  as  a 
function  of  time.  We  perform  an  energy  balance  on  an  element  dR 

assuming  the  streaming  (source)  radiation  is  responsible  for  the  air  burn¬ 
out  at  the  edge  of  the  fireball. 

Energy  Balance  in  dR  in  time  dt 


Here  c  is  the  speed  of  light,  Eg  is  the  energy  density  at  the  edge  of  the 
fireball,  and  all  other  symbols  are  as  previously  defined.  Solving  Eq.  (3) 
for  dR/dt  ,  we  find 


—  -  U 

dt  (U/c)  +  4ttR2Q 

where  we  have  defined 

U  =  4ttR2E  c  . 
s 


(4) 

(5) 


The  boundary  condition  on  Eq.  (4)  is  that  the  fireball  has  a  zero  initial 
radius,  i.  e. , 


R(0)  =  0 


(6) 
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We  now  consider  the  buildup  of  thermal  radiation  in  the  burnout 
sphere.  If  we  assume  the  sphere  to  be  optically  thin,  then  we  can  ignore 
spatial  gradients  and  write  the  rate  equation  for  the  volume  as  a  whole. 

Energy  Balance  for  Thermal  Radiation 


(7) 


The  solution  of  Eq.  (7)  is,  making  use  of  Eq. (16), 


|*R3Et 


=  pxcaTg  Jdt,|ffH3(t')e"pxc(t"t’) 


(8) 


Here  the  new  symbols  introduced  are  ET  ,  the  thermal  radiation  density, 
p  ,  the  air  density,  and  a  ,  the  radiation  constant. 

Finally,  we  consider  the  streaming  radiation.  Since  at  the 
burnout  temperature  the  air  opacity  (x)  is  non-zero,  the  burned 
out  air  will  radiate  and  attempt  to  cool.  We  assume  the  streaming 
radiation  will  be  absorbed  by  the  air  as  it  radiates  to  maintain  its  tem¬ 
perature.  We  perform  an  energy  balance  at  some  radius  r  in  the  in¬ 
terior  of  the  fireball. 


Energy  Balance  for  E„  In  dr  in  time  dt 

1  s  - 


(9) 


i  Jt.  i ,i.i  J«  U  k  -  'JU  -  IfPWURigfgf’ 


|  -  ^  '  jjj|  |pp p| — I — »■■■» 


Equation  (9)  gives  the  result 


d(4frr  E  c)  2  dE 
- g_ - 4rrr  T 


(10) 


Now,  to  maintain  the  air  at  the  burnout  temperature,  the  rate  of  absorption 
of  streaming  energy  must  equal  the  net  rate  of  emission  of  thermal  energy. 
Thus  we  have 


dE 

s 

dt 


pxc(aT 


B 


(ID 


and  Eq.  (10)  becomes 

d(4irr2E  c)  2 

- =  -47rr  pxc(aTg  -ET)  .  (12) 

The  boundary  condition  on  Eq.  (12)  is 

4ur2Es(r,t)c  Yf(t)  .  (13) 


Equation  (13)  is  just  the  energy  conservation  condition  for  an  arbitrarily 
small  sphere  surrounding  the  source. 

We  rewrite  Eq.  (12)  as 


jjp  [4irr2Eg(r,t)c]  =  -47rr2S(t) 


where  we  have  defined 


S(t)  =  pxc[aTg  -ET(t)] 


(14) 


(15) 


Now,  r  in  Eq.  (14)  is  a  path  length  variable  and  t  is  the  associated  time 
at  which  the  radiation  under  consideration  is  at  the  path  length  position  r  . 
Hence  in  Eq.  (14) 


t  -  t  t  +  - 
ref  c 


(16) 
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(17) 


where  tref 


where 


is  any  fixed  reference  time, 
ref  act  c 


is  the  actual  time  of  interest. 


Let  us  set 


Thus  Eq.  (14)  reads 


[4irr2Es(r,t  -~)c]  =  -4irr2S(t  -  ~)  , 


(18) 


where  we  have  dropped  the  subscript  on  t  ^ 
from  r=0  to  r=R  yields 


Integral.  >n  of  Eq.  (18) 


4ffR2Es(R,t)c  =  4»r2Eg(r,t  -  ^)clr=0 

-  j  dr4irr2S(t  -  ^~)  ,  (19) 

0 

or,  using  Eq.  (13), 

H 

4irR2E(R,t)c  =  Yf(t  -  5)  .  J*  dr47rr2S(t  -  ^)  .  (20) 

s  c  Q  c 

It  is  E  (R,t)  that  is  denoted  by  simply  E  in  Eq.  (3), 
s  s 

We  summarize  the  three  pertinent  equations  to  be  solved  for  the 
fireball  radius,  thermal  radiation  density,  and  streaming  radiation  den¬ 
sity  as  a  function  of  time.  These  are  [see  Eqs.  (4),  (8),  and  (20)] 


dR  _ 

u 

(21) 

dt  (U/c)  +  4ttR2Q 

pxcaTg 

Jdt’R3(t,)e"pxc(t't')  , 

(22) 

R3 

0 

Yf(t  -  5)  - 

4?rJ  drr2S(t  -  ^~)  , 

(23) 

o 
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where  [see  Eqs.  (5)  and  (15)] 


U  =  4jtR  E  c 
s 


S  s  pxcfaTg  -  Et]  . 


In  genera),  these  equations  must  be  solved  numerically,  with  starting 


conditions 


R(0)  =  0  , 

Et(0)  =  0  , 


S(0)  =  pxcaTg 

U(0)  =  Yf(0)  , 

dR(0)/dt  =  c  . 


The  numerical  solution  continues  until 


U  =  4irR  E  c  =  0  , 

s  ’ 


i?rR3(Q  +  E  ) 

- — -t - -  =  i  .  02: 

Y  J  dt'f(t') 
o 

The  physical  meaning  of  Eq.  (31)  is  that  the  streaming  radiation  is  entirely 
absorbed  before  it  reaches  the  edge  of  the  fireball.  The  physical  meaning 
of  Eq.  (32)  is  that  the  energy  emitted  by  the  source  is  just  sufficient  to 
maintain  the  air  at  a  temperature  Tg  together  with  the  thermal  radiation 
field  which  has  built  up.  We  shall  discuss  the  finite  difference  analogue  of 
Eqs.  (21)  through  (23)  in  a  later  section  in  this  report. 
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A  quantity  we  shall  find  useful  later  on  is  the  effective  temperature 
of  the  thermal  radiation  field,  defined  by  the  relationship 


aT 


eff 


-  E, 


(33) 


4.  BURNOUT  IN  THE  PRESENCE  OF  A  GROUND 


We  now  modify  the  results  of  the  last  section  to  include  the  in¬ 
fluence  of  the  ground.  The  fireball  will  remain  spherical  until  its  radius 
equals  h  ,  the  height  of  burst.  At  later  times,  we  assume  the  fireball 
is  a  truncated  sphere,  truncated  by  the  air -ground  interface. 

At  such  a  time  when  R>h  ,  we  have  the  schematic 


-y 


Let  us  compute  the  volume  of  various  regions  of  this  sphere.  We  have 
for  all  three  regions  the  obvious  result 


VABC  ’  J*r3  '  <34> 

The  volume  of  regions  A  and  B  is  given  by 

VAR  =  <f  dytrx2  =  J  dy?r(R2-y2)  ,  (35) 

-R  -R 
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(36) 


v  _  w/2o3  „2,  h°v 

VAB  "  nkR  +  R  h  ’  T* 
The  volume  of  region  A  alone  is  given  by 


R 


2ttR 


3  1 


^A  "  2^1  ^rr  J  ^ - 5~  J  > 

o  p  -h/R 


V,  =  |itR2(R  +  h) 


(37) 

(38) 


Now,  we  assume  that  the  yield  directed  into  segment  B  is  r.Dt 
absorbed  in  the  air,  but  is  entirely  absorbed  in  the  ground  and  then  en¬ 
tirely  reradiated  back  into  the  fireball,  all  in  zero  time.  This  energy  is 
then  available  to  increase  the  size  of  the  fireball  in  the  air.  We  account 
for  this  by  increasing  the  yield.  We  let 

Yf(t)  —  Yf(t)d(t)  ,  ,  (39) 

where  d(t)  is  the  enhancement  factor  given  by  the  ratio  of  4tt  to  the  solid 
angle  associated  with  Region  A  .  Thus  we  have 


d(t) 


VABC  2R 
VA  R+h  * 


(40) 


Of  course,  Eq.  (40)  holds  only  for  R  £  h  .  For  R  <  h  ,  we  have 
d(t)  =  1  . 

The  burnout  equations  in  the  presence  of  a  ground  then  become 
[see  Eqs.  (21)  through  (32)] 


dR  _  U 

dt  (U/c)  +  4ttR2Q 


-■y^— -  S  dt'V  (t')e"pKc(t_t,) 
VAB  o 


(41) 


(42) 
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IV 

U  =  Yd(t  -  5.)f(t  -  5)  -  4(T  Jdrr2S(t-  5f£)  , 


with  U  and  S  still  defined  by  Eqs.  (24)  and  (25)  and 


4  r»3 
3*R  ’ 


R  <  h 


(jt(|r3  +  R2h  -  ^)  ,  R  *  h  .  i 

4  3 

In  writing  Eq.  (42)  we  have  replaced  j uR  in  Eq.  (22)  with  VAB  ,  the 
actual  fireball  volume.  The  burnout  phase  of  the  fireball  growth  is  over 
when 


U  =  4jtR  E  c  =  0  , 

s  ’ 


vab(Q  +  et) 

t 

Yjdt’f(t’) 


=  1  . 


The  significance  of  Eqs.  (45)  and  (46)  has  been  discussed  previously. 
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5.  DIFFUSION  PHASE  OF  THE  FIREBALL  GROWTH 


Following  the  end  of  burnout,  we  assume  that  further  fireball 
growth  takes  place  by  a  diffusion  process.  At  the  end  of  burnout,  the 
material  and  thermal  radiation  field  are  in  general  not  in  equilibrium, 
and  we  formulate  the  diffusion  process  to  take  this  into  account.  We 
assume  that  the  diffusion  growth  of  the  fireball  is  described  by  one  di¬ 
mensional  spherical  equation,  but  we  shall  account  for  the  intersection 
of  the  fireball  with  the  ground. 


In  spherical  geometry,  the  energy  conservation  equation  is 


_fi 

St 


(c  T  +  E 
v 


T 


+  E  )  +  4r  (r2F) 
s  2  5r 
r 


S^»(D 

4?rr 


(47) 


where  cv  is  the  air  heat  capacity  (assumed  constant),  F  is  the  radiative 
flux,  and  6(r)  is  the  Dirac  delta  function  indicating  the  source  of  radia¬ 
tion  is  at  the  center  of  the  sphere.  To  obtain  diffusion  theory,  we  assume 
the  radiative  flux  is  proportional  to  the  gradient  of  the  thermal  radiation 
energy  density;  in  particular  we  write 


cX(T)  &ET 
3  fir 


(48) 


where  A(T)  is  the  air  mean  free  path  which,  as  indicated,  depends  upon 
temperature.  We  further  assume  that  A(T)  has  a  cubic  temperature 
dependence,  i.  e. , 


A(T)  =  A(=±-) 
B 


(49) 


where  X  ,  a  constant,  is  the  air  mean  free  path  at  Tg  ,  the  burnout 
temperature.  Use  of  Eqs.  (48)  and  (49)  in  Eq.  (47)  yields 
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(c  T+E-+E  )  = 
5t  v  T  s 


-  2  CX  / 

T  \3  SET 

3  \T 

b)  Sr- 

Yf(t)6(r) 


Finally,  we  assume 

Et  =  eaT4  ,  I 

where  e  is  allowed  to  depend  upon  time  but  not  space.  [If  c  were  also 
allowed  to  depend  upon  space,  Eq.  (51)  could  be  considered  a  definition 
for  e  .  ]  Eq.  (51)  allows  us  to  write 

t3  TF  '  7  £ (af  t7)  =  7  £  ^T*  > 

and  use  of  Eq.  (52)  in  Eq.  (50)  gives  the  diffusion  equation 

4c  <cvT+VEs> -  ?  4c  [f2  £  (t3et>]  +  y  • 

B 


To  proceed  further,  we  form  the  first  two  spatial  moments  of 
Eq.  (53).  Integration  over  all  volume  yields 


O 

A  Jdr47ir  (cvT+Et+Es)  =  Yf(t)  .  (! 

o 

We  now  assume  T,  ET,  and  E  to  be  space  independent  within  the  fire- 

1  s 

ball  and  zero  elsewhere.  This  gives 


If  [  ^R3(cvT+Et+Es)]  =  Yf(t)  . 


4  _,„3 


We  recognize  jffR  as  the  volume  of  a  sphere,  and  interpret  this  volume 
as  the  fireball  volume,  VAB  [see  Eq.  (44)]. 

Thus  Eq.  (55)  becomes 

a  [VAB(cvT+ET+Es)1  =  Y,(t)  •  <5' 
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Equation  (56)  is,  of  course,  just  the  conservation  equation.  Similarly, 

3 

multiplying  Eq.  (53)  by  r  and  integrating  over  all  space,  we  obtain 

±  jdrr3(c  T+E  +E  )  =  _M_jdrrT3E  ,  (57) 

o  V  1  S  21T“L  o  1 
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where  the  right  hand  side  of  Eq.  (57)  resulted  from  two  integrations  by 

parts.  Again  assuming  T,  ET,  and  E  constant  in  space  within  the 

X  s 

fireball  and  zero  elsewhere  yields 

^  [R\T+ VEs>J  =  rV‘et  •  <58> 

For  a  sphere  we  have 


R  = 


(59) 


and  if  we  interpret  the  volume  in  Eq.  (59)  as  ,  the  actual  fireball 

volume,  Eq.  (58)  becomes 


(60) 


Equations  (56)  and  (60)  are  the  first  two  spatial  moments  of  the  diffusion 
equation. 

The  interaction  between  the  material  and  radiation  is  handled  by 
assuming  T,  E^,  and  Eg  satisfy  equilibration  equations  of  the  type 

/T 

i(VABCvT)  ■  M-t)  VAB(ET+VaT  >  ’  (61) 
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(62) 


bt  (VABET)  "  pxc(  T  )  VAB^aT  "ET^ 

3 

i<VABEs>  ■  -P*(t)  VABEs  +  Yf(t)  '  (63) 

The  sum  of  Eqs.  (61)  through  (63)  is  just  Eq.  (56),  the  conservation 
equation.  We  consider  Eqs.  (60)  through  (63)  as  four  equations  for  the 
four  unknowns  V^g  ,  T  ,  E^,  ,  and  Eg  .  These  equations  are  solved 
numerically,  as  we  shall  discuss  in  a  later  section  of  this  report,  with  the 
starting  conditions  on  the  four  dependent  variables  corresponding  to  the 
conditions  at  the  end  of  burnout.  As  before,  it  will  be  useful  to  define  an 
effective  temperature  for  the  thermal  radiation  field  by  the  relationship 


a7 


eff 


=  E, 


(64) 
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6.  FLUX  IMPINGING  UPON  THE  GROUND 


Once  the  fireball  behavior  as  a  function  of  time  has  been  deter¬ 
mined,  it  is  a  simple  task  to  compute  the  flux  of  energy  impinging  upon 
the  ground  at  any  radial  distance  r  (r=0  lies  directly  under  the  device). 

Define  t*  as  the  time  when  the  fireball  first  touches  the  ground 
at  r  .  For  t  <  t*  ,  the  flux  incident  on  the  ground  at  position  r  is 
identically  zero.  The  streaming  radiation  energy  density  at  r  for  t  >  t* 
is  given  by  the  solution  of  the  equation  [see  Eq.  (20)] 


R 


47rR"E  c  =  Yf(t  -  ?)  -  /  d£4rr£2  S(t  -  ^)  , 

s  c  °  ^ 

o 


(65) 


where  R  is  the  fireball  radius  at  t-t*  and  given  by 


-JS  -  2 


R  =  vr  +  h 

and  S  is  given  by  [see  Eq.  (15)] 


(66) 


S  =  pxc(aT^  -  Et) 


(67) 


Then,  the  streaming  flux  per  unit  area  perpendicular  to  the  ground  at  point 
r  is  given  by 


In  addition  to  this  streaming  flux,  the  fireball  thermal  radiation 
field  leads  to  a  flux  in  the  amount 

4 


F  =  cE  h/R 
s  s 


(68) 


FT  aTeff 


(69) 


where  a  =  ca/4  and  Tgf^  is  the  effective  temperature  of  the  fireball 


thermal  radiation  [see  Eqs.  (33)  and  (64)].  Thus  the  total  flux  F 
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impinging  upon  the  ground  as  a  function  of  time  at  some  position  r  is 
given  by 


7.  NUMERICAL  ANALYSIS  OF  FIREBALL 
GROWTH  EQUATIONS 


We  discuss  briefly  the  finite  difference  equations  associated 
with  the  burnout  equations,  (41)  through  (43),  and  the  diffusion  equations, 
(60)  through  (63).  The  finite  difference  equations  conserve  energy  exactly 
and  are  unconditionally  stable  in  the  sense  that  accur  acy  alone  controls 
the  size  of  the  time  step. 


index, 


Equations  (41)  through  (43)  become,  with  i  denoting  the  time 


R.  -  =  R.  +  ra.  (t.  .  -  t.)  , 

l+l  i  'dt'i  i+I  r  ’ 


P’tCaTB  A  1 


JTi  V 


ABi  * 


ZK.-vt 


V  e'pxc(ti+l 'V 

ABj  3 


+  VAB,j+le 


pxc(tifi  -  V]  , 


R  R 

U.  .  =  Yf(t.  ,  -  — )  d(t.  ,  -  -iil) 
i+I  i+I  c  v  i+I  c 


j  =  l 


+  R2  .S 


R.  i  -  R 
i+I 


j+1  \i+l 


[»M« 

Ml  • 


R.  .  -  R.\ 
i+l  3) 

c  / 


i+1  (Uj+1/c)  +  4iR“+1Q 
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I 


. 


3- 

A.(E^  j  +  E  .  ,  —  aT.T.  i) 
1  T,  1+1  s,i+l  1  l+r 


(84) 


.  A.(aA+1  -  Et  .+1) 


(85) 


E±lhmrd  -  -AJ 


t. 


i+1 


iEs,i+l  +  JMi  {  dt’f(t’)  . 


(86) 


Adding  these  equations,  we  find 


t. 


i+1 


(cyT+ET+Es)  =  (cyT+ET+Es)  +  Y  J  dt'f(t')  , 
i+1 


(87 


l  t. 

l 


which  shows  that  the  finite  difference  equations  conserve  energy.  The 
integral  over  f(t)  in  Eqs.  (86)  and  (87)  is  performed  using  a  two-point 
exponential  quadrature  since  f(t)  varies  primarily  exponentially  [see 
Eq.  (1)].  Thus 


J*  dt’f(t’) 

t. 

i 


(At).[f(t.+1)  -  f(t.)] 
^f(ti+1)/f(t.)J- 


(88 


This  same  quadrature  formula  was  used  in  computing  N  ,  the  normali¬ 
zation  integral  in  Eq.  (1),  and  hence  the  integration  of  the  time  dependent 
yield  is  handled  consistently  throughout,  assuring  numeric  energy  con¬ 
servation.  Rearranging  Eqs.  (84)  through  (86),  we  have 


E 


s,  i+1 


t.  i 
l+l 

e0  *  +  Y  r  dt’f(t') 

t. 

'U'K.(tt). - 

1  1 


S,1 


(89 
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(90) 


|t.  +  p 

l  1  cv 

-  ETi 

E  s,  i+l  +  l+A^(At)^ 

(At). 

|  [l+AjUt).] 

1  +  Ai| 

I 

y  c  , 
i  v ) 

|(At)i 

E 


T,  i+1 


Em.  +  A.aT3T.  .(At). 
Ti  1  1  i+l  i 

1  +  A^(At)^ 


Equations  (83)  and  (89)  through  (91)  are  solved  for  the  four  unknowns 
VAB  ■  T  ’  ET  ’  and  Es  ' 


8.  DIFFERENCES  OF  FIREBALL  MODEL  FROM 
EARLIER  (PYATT)  MODEL 


We  list  here  the  major  differences  between  the  fireball  model 
described  in  this  report  and  the  earlier  one  due  to  Pyatt,  et  al. 

a)  The  new  model  is  formulated  for  arbitrary  time  dependence 
of  the  device  output.  A  reasonable  representation  of  this 
time  dependence  is  given  by  Eq.  (1). 

b)  During  burnout,  the  buildup  of  thermal  radiation  is  calcu¬ 
lated  assuming  the  fireball  to  be  optically  thin  [see  Eq.  (7)]. 
The  earlier  model  used  an  optically  thick  equation. 

c)  The  yield  enhancement  factor  [see  Eq.  (40)]  is  more 
physically  reasonable  than  that  used  in  the  Pyatt  model. 

d)  The  earlier  model  uses  an  arbitrary  power  law  to  describe 
the  time  dependence  of  the  fireball  temperature  between  the 
time  the  device  ceases  to  radiate  and  the  time  burnout  is 
over.  No  such  arbitrary  fit  is  used  in  the  current  model. 

e)  The  earlier  model  assumed  equilibrium  between  the  matter 
and  the  thermal  radiation  field  during  the  diffusion  phase  of 
the  fireball  growth.  The  present  model  makes  no  such  as¬ 
sumption,  but  allows  equilibration  to  take  place  simultane¬ 
ously  with  the  diffusion  [see  Eqs.  (61)  through  (63)]. 

f)  The  spatial  dependence  of  the  streaming  radiation  is  cal¬ 
culated  in  the  present  model,  rather  than  assuming  a  l/r^ 
dependence  as  in  the  Pyatt  model  [see  Eq.  (20)]. 

g)  An  ad  hoc  similarity  solution  is  not  used  to  solve  the  diffusion 
equation.  Rather,  a  solution  method  which  correctly  accounts 
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for  both  the  internal  and  radiative  energy  together  is  used 
[see  Eq.  (53)]. 

A  correction  for  the  finite  size  of  the  fireball,  used  in  the 
earlier  model,  is  not  required.  The  thermal  radiation 
field  is  calculated  consistently  for  all  times. 


9.  FINAL  REMARKS 


If  the  ID -2D  synthesis  technique  is  to  be  used  for  future  energy 
coupling  calculations,  it  would  seem  desirable  to  use  the  fireball  model 
discussed  in  this  report.  The  primary  reason  for  this  is  that  this  model 
incorporates  a  more  reasonable  device  output  time  dependence,  and  ana¬ 
lytic  calculations  indicate  that  this  time  dependence  can  significantly 
affect  the  energy  coupled. 

We  intend  to  couple  this  air  fireball  model  with  a  diffusion 
treatment  of  radiative  transfer  in  the  ground,  and  use  this  model  to  per¬ 
form  parameter  surveys  in  the  energy  coupling  area.  This  will  give  us 
a  capability  of  investigating  very  inexpensively  the  sensitivity  of  the  energy 
coupled  to  such  parameters  as  yield,  height  of  burst,  time  dependence  of 
the  device  output,  and  ground  properties.  In  addition,  we  will  be  able  to 
investigate  certain  aspects  of  the  accuracy  of  the  ID -2D  synthesis  tech¬ 
nique,  such  as  the  importance  of  having  a  correct  description  of  the  air 
burnout,  and  the  number  of  radial  slices  required  in  the  ground  for  good 
accuracy. 
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